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基于 Sentinel-2 的 青铜 峡 兴 区 水 稳 和 玉米 种 植 
分 布 早期 识别 
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摘 要 : 及 时 准确 地 掌握 灌区 内 作物 种 植 分 布 对 于 灌溉 水 资源 高 效 配 置 、 农 田 精 准 管理 具有 重要 
导 意 义 。 以 宁夏 青铜 峡 灌区 为 研究 对 象 ,利用 多 时 相 Sentinel-2 卫 星 数 据 ,通过 水 稻 和 玉米 早期 
AE 4 AT ,提取 关键 的 “水 淹 " 信 号 和 “植被 ”信号 LE 
植被 指数 (NDVI) 特 征 值 数 据 集 , 并 通过 样本 分 析 关 键 特征 靖 值 ,构建 水 稻 和 玉米 早期 种 植 分 布 
策 树 模型 ,提取 2022 年 宁夏 青铜 峡 灌 区 水 稻 和 玉米 种 植 的 空间 分 布 。 结 果 表 明 :(1) 玉米 和 水 稳 
期 的 后 半 段 5 月 15 一 31 日 ,水 渡 信 号 和 植被 信号 是 区 分 二 者 关键 时 期 。(2) 基于 早期 作物 物候 特 


EE 时 序 归 一 化 差异 水 体 指数 (MNDWI) 和 和 归 一 


At my 


A H YF c 


的 方法 ,在 5 月 16 一 31 日 获取 的 水 稻 和 玉米 图 像 制图 精度 高 于 90% ,用 户 精度 超过 91% ,总 体 精 度 
超过 90% ,Kappa 系 数 高 于 0.88, 明 显 高 于 同时 期 随机 森林 方法 的 分 类 精度 。(3) 本 研究 提出 的 方法 
在 早期 水 稳 和 玉米 种 植 分 布 提取 方面 具有 较 强 的 适用 性 ,并 且 能 够 在 时 空 尺 度 上 以 较 少 的 实地 样本 
进行 延展 ,同时 在 时 间 上 也 更 有 优势 。 因 此 ,该 方法 为 青铜 峡 灌 区 水 稳 和 玉米 种 植 分 布 早期 调查 提 


供 了 重要 的 方法 支撑 。 
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卫星 遥感 技术 已 被 广泛 地 应 用 作物 识别 分 类 、 
长 势 监测 及 产量 预 佑 等 农业 生产 领域 。 准 确 及 时 
地 获取 作物 种 植 结构 分 布 对 于 农业 水 资源 管理 77 
量 评 佑 农作物 保险 .土地 租赁 供应 链 物流 .灾害 预 
警 以 及 政策 制定 等 方面 具有 极其 重要 的 意义 “"。 以 
往 大 多 数 的 农作物 种 植 分 布 提 取 研 究 ,通常 以 全 年 
多 时 相 琐 感 数据 为 基础 或 作物 全 生命 周期 物候 特 
征 进行 提取 ,只 能 在 年 末 或 第 二 年 度 才能 得 到 作物 
种 植 结 构 分 布 提取 结果 ,时 效 性 较 差 ,无 法 满足 实 
际 生 产 要 求 ”"。 利 用 实时 或 近 实 时 的 卫星 数据 ,在 
夏 灌 之 前 获得 作物 分 类 几 ,可 以 辅助 决策 灌溉 用 水 
量 的 分 配 ,提高 灌溉 的 用 水 效率 * ,减少 因 水 资源 
短缺 之 来 的 产量 损失 风险 。 


收 稿 日 期 : 2023-10-01; 修订 日 期 : 2023-12-28 


基金 项 目 : 国家 自然 科学 基金 项 目 (52209059) ; 宁夏 自然 科学 基金 优秀 青年 项 目 (2023AAC05013 ); 宁夏 回族 自治 区 重点 研发 计划 项 目 


(2021BBF02027) 资 助 


作者 简介 : R46 (1980-7) , 男 ,博士 ,教授 ,主要 从 事 数字 农业 和 智慧 
通讯 作者 : 丁 一 民 (1992-) , 男 ,博士 ,副教授 ,主要 从 事 数 字 农 业 和 


目前 ,已 有 不 少 学 者 关注 到 了 作物 的 早期 识 
别 。 其 中 ,Cai 等 "使 用 Landsat 卫星 2000 一 2015 年 
6 个 光谱 波段 1322 个 多 时 相 场 景 数 据 , 结 合 深度 神 
经 网 络 模 型 进行 可 扩展 的 分 类 过 程 计 算 ;结果 表 
明 ,该 方法 在 同年 7 月 底 实 现 了 对 玉米 和 大 豆 的 分 
类 达到 95% 的 总 体 精度 ,但 该 方法 只 针对 特定 的 农 
作物 进行 分 类 提取 ,在 种 植 结 构 复 杂 的 区 域 还 有 待 
检验 。 郝 鹏 宇 等 利用 2006 一 2013 年 的 MODIS 增 
强 植 被 指数 (EVI) 时 间 序 列 数据 和 农田 数据 ,采用 
免疫 系统 网 络 方法 建立 了 首 苦 、 玉 米 、 高 梁 和 冬 小 
麦 的 参考 EVI 时间 序 列 和 “潜在 ”训练 样本 ;人 研究 结 
果 表 明 , 在 2014 年 作物 分 类 提取 总 体 精度 达到 
94.02% ,接近 使 用 完整 时 间 序 列 数据 的 分 类 结 
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Zhou 5! AJH 4—10 H Foz 250 m MODIS 10 d 合 成 
数据 ,通过 较 短 时 间 序 列 数据 和 一 些 最 佳 变量 组 
合 , 在 加 拿 大 了 萨 斯 喀 彻 温 省 南部 农业 区 进行 作物 种 
植 结构 分 布 提 取 , 但 由 于 MODIS 数 据 空间 分 辨 率 较 
低 ,导致 该 方法 在 田 块 破碎 区 域 对 分 类 精度 造成 较 
影响 。Villa 等 通过 整合 2013 年 和 2014 年 4 一 7 
H 30 m 空 间 分 辩 率 的 Landsat 8 和 X 波 段 SAR 数 据 ， 
基于 分 类 树 方案 ,使 用 多 光谱 数据 集 和 X 波 段 SAR 
数据 派生 的 一 组 季节 性 特征 ,通过 3 种 光谱 指数 和 
X 波 段 SAR 的 植被 和 土壤 条 件 代 理 计 算 来 表达 这 些 
特征 ,最 终 识别 到 生长 季 早 期 7 种 作物 类 型 ,并 通过 
对 2013 年 和 2014 年 2 个 生长 季节 的 验证 ,获得 了 超 
过 87% 的 总 体 精 度 。Azar 5 8H] 2013 年 采集 的 
Landsat 8 0LI 多 时 相 数 据 , 对 意大利 北部 的 7 种 作 
物种 植 结构 进行 分 类 ,该 研究 评 佑 使 用 4 种 监督 算 
法 ,利用 多 时 相 光 谱 指 数 , 如 EVI、 归 一 化 退化 指数 
(CNDFT) ,实现 整体 精度 在 7 月 中 旬 超 过 85% ,在 生长 
季节 结束 的 前 5 个 月 达到 最 高 精度 。 
青铜 峡 灌区 种 植 结 构 复 杂 ,在 生长 季节 同时 存 
在 水 稻 小麦 .玉米 .枸杞 .葡萄 等 作物 ,其 中 玉米 的 
物候 期 和 水 各 高 度 相 似 。 水 淹 信 和 号 是 早期 稻田 的 
一 个 重要 特征 ,但 由 于 漫灌 的 原因 ,玉米 在 早期 也 
存在 偶尔 的 水 淹 人 信号。 尽管 早 期 玉米 和 水 稻 都 存 
在 水 淹 信 号 和 植被 信 号 ,但 是 水 稻 的 “水 淹 - 生 长 ” 
寺 征 更 明显 。 因 此 ,本 研究 利用 多 时 相 Sentinel-2 
数据 来 提高 早期 识别 精度 ,确定 早期 水 称 和 玉米 的 
关键 物候 阶段 和 分 类 特征 ,进而 基于 早期 物候 特征 
方法 提取 水 称 和 玉米 种 植 空间 分 布 ,并 评估 该 方法 
的 潜力 ,为 该 区 域 早期 农 情 调查 和 管理 提供 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 
宁夏 青铜 峡 灌区 (37°30’~38°20’'N, 106?30' - 


106"45' 工 ) 位 于 宁夏 中 北部 , 南 起 青铜 峡 水 利 枢纽 ， 
北 至 石嘴山 y PIER 25 1 ZEE , 东 接 鄂尔多斯 台地 
PUR ,东西 宽 约 50 km, 南 北 长 约 170 km, ERA 
河 冲积 平原 ,包括 青铜 峡 吴忠. 灵 武 . 永 宁 、 银 川 市 
辖 三 区 25 OEY 石嘴山 8 个 市 .县 ,总 面积 约 为 
10245.76 km 。 青 铜 峡 灌 区 地 势 南 高 北 低 ,黄河 流 
经 其 中 ,水 源 充沛 ,土壤 肥沃 ,年 平均 降雨 量 约 
200 mm, 日 照 充 分 , 适 于 发 展 灌溉 农业 ”。20 世 纪 
60 年 代 中 期 建成 了 青铜 峡 水 利 枢纽 工程 ,结束 了 长 
期 以 来 无 坝 引 水 的 历史 ,灌区 设计 灌溉 面积 为 38.8x 
10 hm。 农作物 以 水 稻 、 小麦 .玉米 为 主 , 兼 种 胡 
麻 、 甜 菜 等 ,为 中 国 西北 地 区 重要 产 粮 区 ,被 誉 为 
“塞外 江南 ”。 青 铜 峡 灌 区 主要 农作物 生育 情况 如 
AIR 

1.2 数据 来 源 及 预 处 理 

1.2.1 地 面 验证 数据 为 确保 分 类 和 精度 验证 有 足 
够 数量 和 质量 的 样本 点 ,本 研究 结合 了 野外 实地 调 
查 与 基于 谷歌 地 球 影 像 目 视 解 译 方法 , 共 获 得 了 青 
铜 峡 灌区 内 2022 年 2000 个 有 效 作物 种 植 分 布 样本 
数据 ,其 中 水 稳 641 个 .玉米 616 个 小麦 254 个 其 
他 作物 489 个 。 按 照 3:1 比例 将 所 有 作物 样本 随机 
地 分 为 小 麦 水 稳 、 其 他 作物 (和 葡萄、 枸杞 )、 玉 米 的 样 
本 点 和 分 类 精度 验证 点 (图 1b)。 其 中 ,土地 利用 数 
据 (图 1a) 来 源 于 美国 环境 系统 研究 所 公司 (https:// 
livingatlas.arcgis.com/landcoverexplorer )'', 

1.2.2 Sentinel-2 卫星 数据 及 处 理 Sentinel 系列 卫 
星 是 欧洲 空间 局 (ESA ) 推 出 的 一 组 卫星 , 旨 在 为 地 
球 观 测 提供 开放 高 质量 的 数据 "9 。Sentinel-2 由 2 
颗 卫 星 组 成 , 重 访 周 期 为 Sd,Sentinel-2 卫星 搭载 了 
多 光谱 成 像 仪 。 该 卫星 共有 13 个 波段 ,其 中 4 个 波 
段 的 分 辨 率 为 10 m,3 个 波段 的 分 辨 紊 为 60 m, 光 谱 
覆盖 范围 从 可 见 光 到 近 红 外 至 短波 红外 。 本 研究 
选择 已 经 过 大 气 校正 人 处理 的 Sentinel-2 L2A 级 数 


表 1 青铜 峡 灌 区 主要 农作物 生育 期 情况 汇总 


Tab.1 Summary of main crop growth periods in Qingtongxia irrigation area 
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(a) 土地 利用 类 型 分 布 


>z 


图 例 
号 研究 区 域 
土地 利用 类 型 
mm 水 域 


0 


(b) 野外 调查 样本 点 分 布 
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图 1 青铜 峡 灌 区 土地 利用 类 型 .野外 调查 样本 点 分 布 


Fig. 1 Distributions of land use types and field survey sample points in Qingtongxia irrigation area 


表 2 本 研究 涉及 的 青铜 峡 灌 区 Sentinel-2 
卫星 L2A 级 影像 
Tab.2 Sentinel-2 satellite L2A level images involved in 
this study in the Qingtongxia irrigation area 
获取 时 间 
(年 -月 -日 ) 
播种 前 2022-01-04 ^ S2A-MSIL2A 6 
2022-03-15 S2A-MSIL2A 6 
播种 期 2022-04-06 ^ S2B-MSII2A 5 
2022-04-24 S2A-MSIL2A 5 
2022-05-01 S2A-MSIL2A 6 «8 
2022-05-09 S2B-MSIL2A 6 «5 
HH) 2022-05-16 S2B-MSIL2A 5 «8 
6 
6 
6 
6 
6 
6 


传感器 与 数据 级 别 数量 / 景 ze 


<8 
«5 
«5 
«8 


2022-05-24 S2A-MSIL2A «8 
2022-05-31 S2A-MSIL2A «5 
拔节 期 2022-06-05 S2B-MSIL2A «5 
2022-06-13 S2A-MSIL2A «5 
2022-06-15 «5 
2022-06-30 «5 


S2B-MSIL2A 
S2A-MSIL2A 


据 , 可 以 直接 运用 到 本 研究 中 ""”"。 本 研究 采用 的 
Sentinel-2 卫星 L2A 级 影像 详细 信息 如 表 2 所 示 , 75 
景 影像 中 ,有 74 景 影像 的 云 覆 盖 率 小 于 8%。 其 中 ， 
2023 年 5 月 16 日 的 影像 云 覆 盖 率 23% ,经 过 裁剪 后 


研究 区 部 分 均 为 晴空 状态 ,不 影响 相关 指数 的 计算 。 
1.3 研究 方法 

对 于 玉米 和 水 稻 而 言 ,两 者 早期 的 物候 期 特征 
和 长 势 较为 接近 , 仅 采 用 植被 特征 难以 对 其 进行 
期 区 分 。 因 此 ,本 研究 进一步 融合 水 淹 信 号 开展 玉 
米 和 水 稻 的 区 分 。 然 而 ,早期 玉米 也 存在 水 渡 信 和 号 ， 
特别 是 在 灌溉 后 两 者 水 淹 信 号 也 高 度 相 似 。 但 考虑 
到 玉米 灌溉 后 田间 不 会 持续 有 水 ,本 研究 提出 了 基 
于 多 时 相 早 期 物候 特征 的 水 稻 、 玉 米 种 植 分 布 提取 
方法 。 该 方法 首先 通过 分 析 水 稻 玉米、 小麦 .其 他 
作物 的 年 内 早期 光谱 动态 变化 曲线 ,构建 时 序 归 一 
化 差异 水 体 指数 (MNDWI) 和 归 一 化 植被 指数 (ND- 
VD ,获取 早期 稻田 和 玉米 的 关键 生长 时 期 及 光谱 特 
征 ,并 基于 关键 生长 时 期 光谱 特征 构建 决 分 类 策 树 
模型 提取 水 稻 和 玉米 ,具体 的 技术 流程 如 图 2 所 示 。 
1.3.1. NDVI 和 MNDWI 时 间 序 列 曲线 特征 分 析 水 
稳 在 生长 初期 常 处 于 被 水 淹没 的 生长 状态 ,是 水 域 
与 绿色 植被 的 混合 区 域 ,使 其 更 易于 识别 *“”。 但 
玉米 及 其 他 作物 也 可 能 因 漫 灌 而 存在 水 淹 特 性 ， 
此 本 研究 综合 采用 多 期 遥感 影像 的 MNDWINDVI 
分 别 代表 水 渡 信 号 和 植被 信号 ,对 玉米 和 水 稻 进 行 
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注 :NDVI 为 归 一 化 植被 指数 ;MNDWI 为 归 一 化 差异 水 体 指数 。 下 同 。 


图 2 早期 水 称 \ 玉 米 提 取 技 术 流 程 图 


Fig.2 Flow chart of early rice and corn extraction technology 


早期 识别 。 首 先 , 利 用 NDVI 指 数 将 生育 期 与 水 稳 
和 玉米 有 明显 差异 的 小 麦 .枸杞 及 其 他 作物 进行 易 
除 。 随 后 ,逐步 增加 遥感 影像 数据 ,使 用 MNDWI 和 
NDVI 指 数 进 行 6 次 提取 玉米 和 水 稻 的 分 布 。 
MNDWI 有 效 解决 了 水 域 存 在 阴影 的 问题 , 尤 
其 适用 于 水 体 信 息 的 提取 。 基 于 样本 绘制 的 青 
铜 峡 灌区 主要 作物 的 MNDWI 时 序曲 线 如 图 3 所 示 ， 
具体 的 计算 公式 如 下 : 
MNDWI =(P- P,)/(P,s+P,,) (1) 
式 中 : 已 ,和 已 ,分 别 为 Sentinel-2 卫 星 的 绿 波 段 和 
短波 红外 波段 反射 率 。 
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3 2022 年 青铜 峡 灌 区 主要 农作物 MNDWI、NDVI 时 序曲 线 
Fig.3 Temporal curves of MNDWI and NDVI for major crops in the Qingtongxia irrigation area in 2022 


NDVI 其 数值 变化 与 作物 生长 变化 过 程 相对 
应 。 基 于 样本 绘制 的 青铜 峡 灌区 主要 作物 的 ND- 
VI 时 序曲 线 如 图 3 所 示 , 具 体 的 计算 公式 如 下 : 
NDVI=(P —P,)/(P,s + Pra) (2) 
XP: P P, 分 别 为 Sentinel-2 的 近 红 外 波段 和 
红 波段 反射 率 。 
基于 野外 调查 采集 的 样本 数据 ,计算 青铜 峡 洪 
区 主要 作物 样本 点 在 13 个 时 间 点 ( 表 1、 表 2) 的 ND- 
VI 和 MNDWI 平 均值 3, 构造 青铜 峡 灌 区 主要 农 作 
物 NDVI 和 MNDWI 的 时 间 序 列 曲线 (图 3)。 由 于 水 
稻 和 玉米 的 播种 时 间 相 对 其 他 作物 较 晚 ,其 NDVI 


-0.1 
0 20 40 60 80 100 120 140 160 180 200 
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在 一 开始 明显 低 于 其 他 作物 , 随 着 水 稻 和 玉米 逐渐 
进入 出 苗 期 ,其 NDVI 开 始 逐 步 上 升 。 这 一 现象 构 
成 了 利用 NDVI 来 识别 和 提取 水 稳 和 玉米 信息 的 基 
础 。 水 称 的 MNDWI 自 5 月 上 名 开始 持续 上 升 ,直至 
5 月 下 旬 达 到 峰值 ,该 峰值 稳定 在 0.08-0.12 之 间 。 
值得 注意 的 是 ,在 水 稻 生 长 初期 ,其 MNDWI 的 变化 
幅度 显著 高 于 玉米 .小 麦 以 及 其 他 经 济 作物 ,如 枸 
杞 .葡萄 等 。 这 一 显著 差异 为 利用 MNDWI 有 效 区 
分 水 稻 和 其 他 作物 ,尤其 是 玉米 ,提供 了 重要 基础 。 
13.2 基于 决策 树 的 玉米 和 水 稳 早 期 识别 模型 基 
青铜 峡 灌区 的 种 植 特点 ,本 研究 构建 了 一 个 可 动 
态 更 新 的 决策 树 模型 ,以 遥感 影像 为 数据 源 , 随 着 
遥感 数据 的 融入 不 断 优 化 识别 精度 ,实现 对 玉米 和 
水 稻 的 早期 精准 识别 。 首 先 ,利用 5 月 1 日 的 NDVI 
图 像 和 土地 利用 数据 , 噜 除了 灌区 内 生育 期 较 早 的 
Ft E Cn BER o INE .枸杞 .葡萄 等 )。 模 型 进 一 
步 利 用 后 续 6 个 关键 时 期 的 作物 植被 信号 和 水 淹 信 
号 特征 ,通过 样本 点 概率 分 布 统计 来 确定 区 分 玉米 
和 水 稻 的 分 类 阔 值 。 有 具体 而 言 ,模型 在 5 月 1 日 设 
定 NDVI<0.106 ,以 进一步 精确 玉米 和 水 稻 的 种 植 区 
域 。 在 此 基础 上 ,使 用 5 月 16 A MNDWI>-0.08 对 
水 稻 和 玉米 进行 初步 区 分 ,得 到 分 类 结果 一 。 随 着 
作物 进入 出 苗 期 ,植被 信号 日 益 明 显 , 采 用 5 月 24 
日 0.125<NDVI<0.675 和 5 月 31 H 0.025<NDVI< 
0.575 以 及 相应 MNDWI 的 分 类 阔 值 MNDWTI>-0.08 
和 MNDWI-0.1 进 行进 一 步 区 分 ,依次 得 到 分 类 结 
果 二 和 三 。 在 作物 生长 的 后 续 阶 段 ,模型 继续 根据 
NDVI 和 MNDWI 的 变化 特征 来 提取 水 稻 和 玉米 的 
种 植 分 布 。 在 分 类 结果 三 的 基础 上 ,通过 设 定 6 月 5 
日 -0.125<NDVI<0.375 fil MNDWIz- 0.08 £5 [E] fii , 
得 到 分 类 结果 四 。 进 入 拔节 期 ,由 于 2 种 作物 的 植 
被 信号 更 加 相似 ,模型 转 而 利用 水 淹 信 号 的 差异 ， 
取 6 月 13 日 MNDWI=-0.15 进行 区 分 ,得 到 分 类 结 
果 五 。 最 后 ,在 6 月 下 旬 , 当 水 稳 和 玉米 的 水 淹 信 号 
差异 相对 较 小 时 ,采用 6 月 30 日 0.125<NDVI<0.675 
fil MNDWIZ-0.3 的 准 值 进行 最 终 区 分 ,得 到 分 类 结 
果 六 。 随 着 作物 生长 ,不 断 通 过 阔 值 和 条 件 判 断 ， 
构建 了 一 个 高 效 的 决策 树 分 类 模型 (图 4)。 
1.3.3 基于 随机 森林 算法 的 作物 分 类 随机 森林 算 
法 最 早 由 Breiman 提出 号 ,是 一 种 有 监督 的 机 器 学 
习 算 法 ,监督 学 习 是 机 器 学 习 的 一 个 子 类 。 这 种 类 


识别 图 像 中 的 不 同类 型 。 随 机 森林 由 于 其 准确 性 ， 
简单 性 和 灵活 性 , 现 已 成 为 最 广泛 使 用 的 分 类 算法 
之 一 六。 随机 森林 通过 集成 学 习 的 方式 ,充分 发 
挥 了 决策 树 的 优势 ,并 通过 随机 特征 选择 和 投票 机 
制 来 提高 模型 的 泛 化 能 力 和 稳定 性 。 由 于 随机 森 
林 属 于 机 器 学 习 工 具 ,其 分 类 算法 主要 利用 了 影像 
的 光谱 人 信息。 因此 ,其 样本 点 因 遵 循 “质量 胜 于 
数量 "原则 , 即 样本 点 要 保证 较 高 的 精度 。 本 研究 
中 随机 森林 算法 输入 的 影像 和 调查 样本 与 上 述 决 
策 树 方法 一 致 ,并 同样 按时 间 顺 序 ,逐步 增加 数据 
进行 了 6 次 识别 (图 5)。 

13.4 水 稻 、 玉 米 提 取 方 法 精度 评估 本 研究 采用 
混淆 矩阵 以 及 计算 制图 精度 、 用 户 精 度 、 总 体 精度 和 
Kappa 系数 对 水 稳 和 玉米 提取 进行 精度 评定 …”。 


2 结果 与 分 析 


2. 水 稳 和 玉米 提取 精度 评价 

本 研究 依托 早期 物候 特征 对 水 稻 和 玉米 的 种 
植 分 布 进行 了 综合 分 类 ,总体 精度 均 高 于 81% 以 
上 。 尤 其 值得 注意 的 是 ,在 出 苗 期 ,对 2 种 作物 的 制 
图 精度 均 超 过 93% ,用 户 精 度 均 超过 91% ,Kappa 系 
数 均 达 到 0.8 以 上 ( 表 3), 具 体 的 分 类 结果 详 见 图 
6。 然 而 , 随 着 作物 的 生长 进入 拔节 期 ,分 类 精度 呈 
下 降 趋势 ,在 拔节 期 后 期 总 体 分 类 精度 降低 至 
81.0671% , Kappa 系数 降 至 0.609。 这 主要 是 由 于 随 
着 作物 的 生长 导致 地 表 被 植被 遮盖 ,使 得 水 淹 信 号 
在 区 分 水 稳 和 玉米 时 逐渐 减弱 ,进而 导致 一 定 程度 
的 误 分 类 现象 。 上 述 结果 表明 ,5 月 16—31 日 即 出 
苗 期 的 植被 信号 和 水 淹 信 号 ,是 区 分 水 稳 和 玉米 的 
最 佳 时 间 窗 口 。 
2.2 不 同 分 类 方法 提取 精度 表现 对 比 

基于 随机 森林 分 类 方法 的 作物 提取 精度 如 表 4 
所 示 ,通过 表 3 和 表 4 精 度 对 比 可 以 发 现 ,基于 早期 
物候 特征 的 总 体 精 度 较 高 ,尤其 是 在 出 苗 期 ,用户 
精度 和 制图 精度 上 表现 的 更 明显 ,用 户 精度 显著 高 
tH 2%~4% ,制图 精度 总 体高 出 4%~7% , 即 基 于 早期 
物候 特征 提取 水 稻 和 玉米 的 错 分 误差 和 漏 分 误差 
更 小 。 基 于 早期 作物 物候 特征 方法 提取 水 稻 和 玉 
米 时 ,能 更 好 地 避免 因 出 现 偶尔 水 淹 信 和 号 而 导致 的 
错 分 。 相 比 之 下 本 研究 提出 的 分 类 方法 ,充分 利用 


型 的 学 习 依赖 于 分 类 标签 ,以 生成 一 个 函数 模型 来 


了 作物 生长 早期 物候 期 特征 , 既 兼 顾 考 虑 了 水 稻 和 
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5 日 .6 月 13 日 .6 月 30 日 的 


31 日 .6 月 5 日 .6 月 13 日 .6 
月 30 日 的 归 一 化 差异 水 体 指数 。 
图 4 分 类 决策 树 规 则 
Fig.4 Rules of classification decision tree 
em | 玉米 早期 的 植被 信号 特征 ,同时 也 考虑 了 两 者 之 间 
水 淹 信 号 的 差异 ,大 大 提高 了 分 类 精度 。 
: 图 7 为 选取 的 2 个 典型 区 域 分 类 结果 ,以 2022 
REM 年 5 月 16 日 (图 7a) 和 9 月 17 日 (图 7b) 的 谷歌 地 球 
TEY A il we 分 
BM E 卫星 影像 数据 为 参考 ,观察 到 2 种 方法 的 分 类 差 


Fig. 5 Flow of random forest algorithm 异 。 在 1 写 区 域 ,随机 森林 方法 误 将 一 部 分 水 体 区 
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表 3 基于 早期 物候 特征 方法 的 精度 评价 


Tab.3 Accuracy evaluation based on early phenological feature extraction methods 


玉米 TAR TER " 

wali RM 制图 精度 ” ”用户 精度 /% RB Erie BEI Kappa 系数 
eA 5 月 1 日 一 16 日 95.62 94.33 93.02 92.66 97.0413 0.920 
出 苗 5 月 1 日 一 24 日 95.25 93.31 94.87 93.34 95.4182 0.908 
5 月 1 日 一 31 日 93.49 92.87 93.54 91.56 91.5336 0.882 
拔节 5 月 1 日 一 6 月 05 日 90.96 91.41 92.06 93.14 90.6971 0.851 
5 月 1 日 一 6 月 13 日 89.36 90.61 90.84 91.33 89.5663 0.801 
5 月 1 日 一 6 月 30 日 63.53 92.86 95.87 75.69 81.0671 0.609 

(a) 5 月 16 日 (b) 5 月 24 日 (c) 5H31H 


图 例 _ | mam EK B Gc 


6 不 同方 案 提 取 的 水 称 玉米 种 植 空间 分 布 


Fig. 6 Spatial distributions of rice and corn cultivation extracted with different schemes 


表 4 基于 随机 森林 分 类 方法 的 精度 评价 


Tab. 4 Accuracy evaluation based on random forest classification method 


玉米 水 稳 ae N 

oak a MARE HRE RIM HPE 总体 精度 /% Kappa t 
JE RIE 5H16H 89.34 87.56 90.01 89.35 90.01 0.89 
ois 5H24H 89.98 90.20 89.97 90.02 90.06 0.87 

5H31H 87.34 88.75 86.37 89.23 89.09 0.84 
拔节 6 H05H 90.02 90.51 90.07 91.24 90.76 0.90 

6H 13H 88.09 90.03 89.67 90.01 89.12 0.85 

6 H 30 H 91.03 90.5 86.08 87.09 86.57 0.86 
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(a) 1 号 区 域 


(b) 2 号 区 域 


(c) 1 号 区 域 随机 森林 的 分 类 结果 


(e) 1 号 区 域 基 于 早期 物候 特征 的 分 类 结果 


( 2 号 区 域 基于 早期 物候 特征 的 分 类 结果 
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m 


KAC 典型 区 域 pM 于 玉米 


图 7 典型 区 域内 2 种 分 类 方法 效果 对 比 


Fig.7 Comparison of the effectiveness of two classification methods in typical regions 


域 错 分 为 水 稻 ( 图 7c)。 在 2 号 区 域 ,随机 森林 方法 
漏 分 了 部 分 玉米 地 (图 7d)。 相 较 于 随机 森林 分 类 
方法 ,基于 本 研究 提出 的 基于 早期 物候 特征 的 分 类 
方法 在 1 号 区 域 和 2 号 区 域 避免 了 错 分 (图 7e) 和 漏 
分 (图 7f) 情 况 。 因 此 ,基于 早期 物候 特征 的 分 类 方 
法 展现 出 了 更 高 的 准确 性 和 可 靠 性 。 
2.3 不 同 分 类 方法 在 早期 识别 中 的 表现 

基于 早期 物候 特征 方法 和 随机 森林 方法 的 Kap- 
pa 系数 最 大 值 分 别 为 0.92(2022 年 5 月 16 日 ) 和 0.90 
(2022 年 6 月 5 日 ), 在 5 月 31 日 之 前 ,基于 早期 物候 
村 征 方 法 的 整体 精度 较 高 ,并 且 Kappa 系 数 首次 大 于 
0.9 的 日 期 更 早 ( 表 5$)。 这 个 时 期 水 稻 和 玉米 的 水 淹 
信号 差异 显著 ,使 得 早期 物候 特征 的 分 类 方法 表现 
出 更 大 的 优势 。 需 要 强调 的 是 ,在 本 研究 的 分 类 方 
法 中 ,总体 精 度 呈 逐渐 降低 的 趋势 ( 表 3)。 根 据 图 3 


#5 Kappa 系数 分 析 
Tab.5 Kappa coefficient analysis 


i j K ARCO. 
i ee e E ane 9 
基于 早期 物候 特 092 2022-05-16 2022-05-16 
征 方法 
随机 和 森林 方法 0.90 2022-06-05 2022-06-05 


所 示 ,在 2 种 作物 分 类 中 ,早期 水 稻 MNDWI 变 化 最 
大 ,尤其 是 在 出 茵 期 差异 较 大 ,因此 在 这 个 时 间 段 内 
分 类 精度 相对 较 高 。 然 而 , 随 着 6 月 初 水 称 进 入 拔节 
期 ,水 淹 信 号 逐渐 降低 ,导致 分 类 精度 逐渐 下 降 , 这 
现象 需要 在 未 来 的 研究 中 进行 深入 分 析 和 改进 。 


3 讨论 


本 研究 以 Sentinel-2 时间 序 列 数据 为 基础 , 采 
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用 早期 物候 特征 方法 评估 青铜 峡 灌区 早期 水 稳 和 
玉米 的 分 类 效果 ,并 与 随机 森林 分 类 方法 进行 对 比 
分 析 。 实 验 验 证 本 人 研究 方法 在 早期 水 稳 和 玉米 提 
取 方 面 的 有 效 性 ,其 分 类 精度 在 5 月 底 达 到 91% , 满 
足 了 在 夏季 灌溉 之 前 获取 作物 分 类 图 的 要 求 。 此 
外 , 相 较 Cai 等 "方案 ,本 研究 使 用 更 少数 据 ,取得 
更 高 分 类 精度 , 且 时 间 更 早 。 在 数据 选择 上 ,本 研 


玉米 提取 方法 。 当 NDVI-0.125 时 ,被 认定 为 水 稳 
和 玉米 的 植被 信号 区 域 。 在 植被 信号 提取 的 区 域 
范围 内 , 当 MNDWI=-0.08 时 ,被 判定 为 水 稻 ; 反 之 ， 
则 为 玉米 。 该 方法 在 5 月 16 一 31 日 期 间 获 取 的 水 
稻 和 玉米 图 像 制图 精度 高 于 90% ,用 户 精度 超过 
9196 ,总 体 精度 超过 90% , Kappa 系数 超过 0.88。 这 
一 方法 有 望 满足 辅助 决策 青铜 峡 灌区 灌溉 的 用 水 


究 充 分 利用 了 Sentinel-2 数 据 的 时 间 和 空间 分 辨 率 
优势 ,尤其 在 田 块 面积 较 小 区 域 ,相对 于 郝 鹏 宇 等 * 
研究 ,该 方法 在 作物 分 类 提取 上 表现 更 为 出 色 。 与 
Villa 等 研究 相 比 ,本 研究 采用 较 少 的 植被 指数 ， 
同时 该 方法 更 具 物 理解 释 性 , 量 在 时 间 方 面 更 具 优 
势 。 

虽然 本 研究 已 得 到 较 高 精度 分 类 结果 ,但 是 依 
然 存 在 以 下 不 足 。 首 先 , 基 于 作物 早期 物候 特征 的 
植被 指数 阔 值 分 类 方法 可 能 导致 分 类 图 斑 破 碎 化 ， 
尤其 在 田 块 分 类 完整 性 上 可 能 显得 不 连续 和 破 
雄 。 未 来 的 研究 可 以 考虑 在 基于 田 块 尺度 上 进行 
作物 分 类 提取 , 以 提高 分 类 结果 的 一 致 性 。 其 次 ， 
本 研究 在 利用 早期 植被 信号 和 水 淹 信 号 进行 水 稻 、 
玉米 种 植 分 布 提取 时 , 随 着 作物 生长 后 期 ,水 淹 信 
号 逐渐 减弱 ,导致 后 期 的 分 类 精度 下 降 。 在 以 后 的 
研究 中 应 考虑 针对 不 同 生 长 阶段 采用 更 为 灵活 的 
特征 进行 分 类 ,以 应 对 作物 生长 过 程 特征 变化 的 复 
杂 性 。 同 时 ,考虑 到 不 同 地 区 的 气候 和 种 植 管理 方 
式 的 不 同 ,研究 者 在 使 用 该 方法 时 需要 结合 本 地 具 
体 作 物 生 长 物候 期 ,确定 早期 提取 水 稻 和 玉米 种 植 
分 布 的 最 佳 时 间 。 最 后 ,本 研究 主要 结合 作物 早期 
的 物候 特征 和 多 时 相 光 谱 特 征 进行 作物 种 植 分 布 
提取 ,但 不 可 避免 地 出 现 了 错 分 和 漏 分 的 现象 。 在 
未 来 的 研究 中 ,可 以 考虑 引入 深度 学 习 方 法 ,充分 
利用 其 在 遥感 图 像 分 类 方面 的 卓越 自主 学 习 和 特 
征 提 取 能 力 。 通 过 结合 深度 学 习 技术 ,利用 影像 的 
其 他 信息 ,如 耕地 边界 .作物 纹理 特征 、 作 物 播种 方 
SUSE OR BUE TRE, EREU RERA 
进一步 发 展 提供 新 的 前 景 。 


4 结论 


(1) 早期 水 稻 和 玉米 提取 的 关键 物候 期 及 分 类 
特征 为 :5 月 16—31 日 的 水 渣 和 植被 信号 为 区 分 水 
稻 和 玉米 种 植 分 布 的 关键 物候 期 , 即 出 苗 期 后 期 。 

(2) 建立 了 基于 早期 关键 物候 期 特征 的 水 稻 和 


效率 ,减少 由 于 水 资源 短缺 导致 的 产量 损失 风险 。 

(3) 相对 于 随机 森林 的 分 类 成 果 ,本 研究 提出 
的 基于 早期 关键 物候 期 特征 的 方法 不 仅 在 分 类 准 
确 度 上 表现 更 为 出 色 ,而 且 在 时 间 上 能 够 更 早 地 获 
取 作 物种 植 空间 分 布 等 方面 展现 出 明显 的 优势 。 
特别 值得 强调 的 是 ,在 北方 单 季 种 植 区 ,该 方法 具 
备 广泛 的 适用 性 ,可 被 广泛 推广 应 用 。 
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Early identification of rice and corn planting distribution in Qingtongxia 


irrigation area based on Sentinel-2 
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of Digital Water Network Management of the Yellow River in Ningxia Hui Autonomous Region, Yinchuan 750021, Ningxia, China) 


Abstract: Timely and accurate understanding of crop distribution within irrigation areas is essential for the effi- 
cient allocation of irrigation water resources and precise field management. This study focuses on the Qingtongx- 
ia irrigation area in Ningxia, China, employing multitemporal Sentinel-2 satellite data to analyze early character- 
istics of rice and maize. Key "flooding" and "vegetation" signals are extracted, and a time-series dataset compris- 
ing the modified normalized difference water index (MNDWI) and normalized vegetation index (NDVI) is con- 
structed. By analyzing sample thresholds for these key features, a decision tree model for the early planting distri- 
bution of rice and maize is established, facilitating the extraction of the spatial distribution for rice and maize 
planting in the Qingtongxia irrigation area in 2022. The results reveal the following: (1) During the latter half of 
the maize and rice seedling stages, from May 15 to 31, flooding and vegetation signals are crucial for differentiat- 
ing between the two crops. (2) Based on the early crop phenological characteristics, the mapping accuracy of rice 
and corn images obtained from May 16 to May 31 was higher than 90%, with user accuracy exceeding 91% and 
overall accuracy exceeding 9096. The Kappa coefficient was higher than 0.88, significantly higher than the classi- 
fication accuracy of the random forest classification method during the same period. (3) The proposed method 
demonstrates strong applicability in the early extraction of rice and maize planting distribution, requiring fewer 
ground samples for extension across both spatial and temporal scales. Therefore, this method provides significant 
support for early investigations of rice and maize planting distribution in the Qingtongxia irrigation area. 

Key words: Qingtongxia irrigation area; Sentinel-2; normalized vegetation index; normalized difference water 


index; decision tree; rice; maize 


